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ABSTRACT 

Recent observations of galaxy clusters in radio and X-ray indicate that cosmic rays and magnetic 
fields may be energetically important in the intracluster medium. According to the estimates based 
on theses observational studies, the combined pressure of these two components of the intracluster 
medium may range between 10% ^ 100% of gas pressure, although their total energy is probably time 
dependent. Hence, these non-thermal components may have influenced the formation and evolution of 
cosmic structures, and may provide unique and vital diagnostic information through various radiations 
emitted via their interactions with surrounding matter and cosmic background photons. We suggest 
that shock waves associated with cosmic structures, along with individual sources such as active galactic 
nuclei and radio galaxies, supply the cosmic rays and magnetic fields to the intracluster medium and to 
surrounding large scale structures. In order to study 1) the properties of cosmic shock waves emerging 
during the large scale structure formation of the universe, and 2) the dynamical influence of cosmic 
rays, which were ejected by AGN-like sources into the intracluster medium, on structure formation, we 
have performed two sets of N-body/hydrodynamic simulations of cosmic structure formation. In this 
contribution, we report the preliminary results of these simulations. 
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I. INTRODUCTION 

Formation and evolution of cosmic structures in¬ 
cluding clusters of galaxies have been driven by the 
gravity of matter, and the matter continues to ac¬ 
crete onto the cosmic structures as a part of hier¬ 
archical structure formation in the cold dark matter 
paradigm. The accreting matter has a typical velocity 
up to ~ a few lO^km s~^, and its baryonic component 
ends as accretions shocks of Mach number from a few to 
a few hundred (see, e.g., Ryu & Kang 1997a). The grav¬ 
itational energy of the gas is dissipated via collisionless 
shocks, heating the intracluster gas to 10® — 10®K. Cos¬ 
mic shocks, along with the hierarchical clustering of 
predominant dark matter, generate complex flow pat¬ 
terns of gas inside cosmic structures, which, in turn, in¬ 
duce weaker, internal shocks. As a result, shock waves 
formed during the large scale structure formation span 
a wide range of scales from a few kpc to a few Mpc, 
and have a wide range of Mach numbers from one to 
a few hundred (see, e.g., Ryu & Kang 1997b, Miniati 
et al. 2000). The shocks can heat up the gas to the 
temperature ranging from ~ 10^ — 10® K in pancakes 
up to ~ 10® K in clusters (see, e.g., Kang et al. 1994, 


Cen & Ostriker 1999). 

In addition to heating gas, the shock waves are 
likely to accelerate high energy particles (cosmic rays) 
through the so-called diffusive shock acceleration, pro¬ 
vided that there exist weak magnetic fields in the in- 
tergalactic space (see, e.g., Blandford & Eichler 1987, 
Kang, Ryu, & Jones 1996). Extended regions popu¬ 
lated by cosmic ray electrons have been observed in 
some clusters for more than thirty years through dif¬ 
fuse, nonthermal radio emissions (see e.g., Kim et al. 
1989). Although cosmic ray protons produce y-rays 
via TT® decay following inelastic collisions with gas nu¬ 
clei, such y-rays have not yet been detected (Sreekumar 
et al. 1996). The observations of nonthermal radia¬ 
tions from cosmic ray electrons in clusters, however, 
suggest that cosmic ray protons could be energetically 
important, and possibly in energy equipartition with 
gas. Cosmic shock waves may be the main sources of 
such cosmic rays, while active galactic nuclei including 
radio galaxies should also provide a significant amount 
of cosmic rays to the intracluster medium. 

Cosmic shock waves could serve also as sites for gen¬ 
eration of weak seeds of magnetic fields by the Bier- 
mann battery mechanism. It was proposed that these 
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seeds could be amplified to strong magnetic fields of up 
to a few /zG in clusters, if flows there can be described 
as the Kolmogoroff turbulence (Kulsrud et al. 1997). 
Although further development into coherent magnetic 
fields is unclear, since there is as yet no detailed theory 
capable of describing this process, observations suggest 
the existence of cluster magnetic fields of a few /xG 
strength and ~ 10 kpc coherent length (see, e.g., Glarke 
et al. 2001). 

We have studied the properties of cosmic shock 
waves and the dynamical roles of cosmic rays in clusters 
by high-resolution cosmological hydrodynamic simula¬ 
tions. For the latter, we adopted an AGN (active galac¬ 
tic nucleus) origin model in which cosmic ray energy 
is deposited into the intracluster medium by AGNs, 
while the cosmic rays accelerated at cosmic shocks are 
ignored. This is because there is yet no practical nu¬ 
merical scheme which can follow the non-linear diffusive 
shock acceleration of cosmic rays in multi-dimensional 
simulations. In the next section we describe a cosmo¬ 
logical hydrodynamic simulation and discuss the sta¬ 
tistical properties and roles of shock waves in cosmic 
structures. The final section briefs simulations which 
include the dynamical influences of cosmic rays in the 
course of formation and evolution of cosmic structures. 

II. SHOCKS FROM A HIGH-RESOLUTION 
COSMOLOGY SIMULATION 

The properties of shock waves associated with cos¬ 
mic structures have been studied through a high-resolu¬ 
tion cosmology simulation. For the simulation, an 
Eulerian hydro-|-N-body cosmology code (Ryu et al. 
1993) has been used. The cold dark matter cosmology 
with a cosmological constant (ACDM) has been em¬ 
ployed with the following parameters: flsM = 0.043, 
^DM = 0.227, and = 0.73 {^bm + = I), 

h = i7o/(100 km/s/Mpc) = 0.7, and erg = 0.8. A cu¬ 
bic region of size 100ft-“^ Mpe at the current epoch has 
been simulated inside a computational box with 1024^ 
gas and gravity cells and 512^ dark matter particles, 
allowing a spatial resolution of 97.7h~^ kpc. The sim¬ 
ulation is the largest of this kind to date. 

Fig. 1 illustrates typical cosmic structures found in 
the simulation, showing bremsstrahlung X-ray emissiv- 
ity (upper-left), gas density (upper-right), gas temper¬ 
ature (lower-left) and shock Mach number (lower-right) 
distributions. It represents a slice of (25/i“^Mpc)^ cen¬ 
tered on a hot cluster in the simulation box. While X- 
ray clusters are scattered, gas density and temperature 
reveal connected structures mostly through pancakes 
and filaments. The outer bounds of shock waves fol¬ 
low closely those of gas temperature distribution. But 
Mach number distribution shows rich, complex network 
of weak shock waves inside the filaments bounded by 
strong accretion shocks. Note that among the protrud¬ 
ing structures around the central cluster, thin features 
in T and M distribution with a thickness of ~ lh~^ 
Mpe correspond to cross sections of pancakes, while 


somewhat thick features with a thickness of ~ 3 — 5/i ^ 
Mpe belongs to filaments. 

Fig. 2 shows three-dimensional perspective of cos¬ 
mic structures through volume rendering of gas tem¬ 
perature. Top panel covers the whole simulation box, 
while lower two panels represent two different perspec¬ 
tives for a (25/i“^Mpc)^ portion centered on the same 
cluster shown in Fig. 1. Images show clearly three 
topological features, i.e., pancakes, filaments and knots 
(clusters) in cosmic structures. Filaments lie on pan¬ 
cakes, and clusters are distributed along filaments or 
positioned at the intersections of filaments. 

Cosmic shocks are most commonly characterized as 
either accretion shocks, if they result from infall of dif¬ 
fuse, intergalactic gas onto the perimeter of clusters, or 
merger shocks, if they result from collisions of two clus¬ 
ters. A quick glance at Figs. 1 and 2 reveal that this is 
an overly simplified picture. The hierarchical cluster¬ 
ing process for structure formation produces extremely 
complex shock structures inside, around and outside 
clusters. These shock waves are neither spherical nor 
identifiable by simple surfaces. Indeed, they intersect 
each other, forming nested shock surfaces, extending 
from outer bounds of pancakes and penetrating deep 
inside clusters. In addition, collisions between flows 
in filaments can lead to shocks, and accretion shocks 
are often hard to be distinguished from merger shocks, 
given the complexity that accompanies the accumula¬ 
tion of mass in regions where clusters are forming. Dis¬ 
sipation at these shocks provides the basic heating of 
the intracluster medium, although other processes, in¬ 
cluding feedback from active galactic nuclei and galax¬ 
ies may also be important contributors. 

According to Miniati et al. (2000) and our analysis 
of the current simulation, cosmic shocks come in wide 
ranges of scales and strengths. Shock waves have been 
identified inside, around and outside clusters on scales 
from ^100 kpc to a few Mpe. Shock Mach number de¬ 
pends on the gas temperature, since flow velocities tend 
generally to be of order 10^ km s“^ both inside and out¬ 
side clusters. Since shocks are found more frequently 
inside clusters, where gas temperature is high (T ^ 10® 
K), there are more shocks with smaller Much num¬ 
ber. The shock surface-area distribution shows a peak 
at M = 1. However, to the extent that the shocked 
gas has been virialized, moderate strength shocks with 
Mach number roughly in the range 2 ^ M ^ 3 would 
contribute most to heating of the intracluster medium. 

Gosmic shocks are capable of accelerating cosmic 
rays efficiently, so deserve a close scrutiny in that re¬ 
gard. It is important to remember that, since cosmic 
rays are effectively tied to the intracluster medium up 
to pretty high energies, the integrated shock history 
of the intracluster medium determines the character of 
the cosmic ray population. So, one would not expect 
cosmic rays produced in this way to be only associated 
with recent merger events, for example. Our analysis 
of the simulation indicates the shocks, which contribute 
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Fig. 2 .— Three-dimensional volume rendering of gas temperature in the ACDM simulation with 1024® grid cells. Top: The 
box of (100/i~®Mpc)®. Middle and Bottom: Two different perspective views of a (25/i~®Mpc)® portion which encompasses 
the slice in Fig. 1. 
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most to acceleration of cosmic rays, would be those in 
the range 3 ^ M ^4. Miniati et al. (2001) and our 
analysis show such shock waves can accelerate cosmic 
ray protons with pressure up to ~ 50% of gas thermal 
pressure. 

Details of the simulation in this section will be re¬ 
ported elsewhere (Ryu et al. 2002b; Kang et al. 2002). 


III. COSMOLOGY SIMULATIONS WITH 
COSMIG RAY DYNAMIGS 


The simulations used to study the dynamical effects 
of cosmic rays have also employed the ACDM cosmol¬ 
ogy with the parameters same as those in the previ¬ 
ous section, and have been performed with the same 
cosmology code. However, a cubic region of smaller 
comoving size, 75/i“^ Mpc, has been chosen for com¬ 
putational domain, and smaller numbers of cells and 
particles, 512^ cells for gas and gravity and 256^ parti¬ 
cles for dark matter, have been used. Two simulations 
have been done, one with cosmic ray dynamics included 
and the other without cosmic rays. 

The following equation for cosmic ray pressure has 
been solved in addition to the usual standard set of 
equations for dark matter and gas (Ryu et al. 1993), 


dpcR , 1 dPcR , 'YcR D 9uk 

—I—-1- ^CR ^— 

ot a oxk CL oxk 


-3(7cfl-l)-pci?- 


Here, a is the expansion parameter, Xk is the comoving 
length, pcR is the comoving pressure of cosmic rays, 
Uk is the proper peculiar velocity, and ^CR is the adia¬ 
batic index of cosmic rays which has been assumed to 
be 4/3. Diffusion of cosmic rays has been ignored, since 
the diffusion length is much shorter than the computa¬ 
tional cell size for most of cosmic rays. The dynamical 
feedback of cosmic rays to gas has been incorporated by 
including the {1/a){dPc r/ dxu) term in the momentum 
equation, along with the {1/a)uk{dPcR/dxk) term in 
the energy equation. 

It has been assumed that the sources, which deposit 
cosmic ray energy into the intracluster medium, form 
at 40 different epochs after the redshift 2 : = 10, if the 
following criteria are satisfied in each grid cell 


Mgas > 


3 X 10^° 
1 + z 


h ^Mq, 


duk 

dxk 


0 , 


( 2 ) 


where Mgas is the total gas mass inside the cell. It has 
been further assumed that each source ejects the follow¬ 
ing amount of cosmic ray energy into the intracluster 
medium 

Ecr = 3 X lO^h-^Mg X c^. (3) 

Note that at z = 0 this translates into the cosmic ray 
efficiency, EcR/MgasC^ = 10“®. 

Bottom panel of Fig. 3 shows that the source for¬ 
mation history, which has been realized by the criteria 
(2), increases first, and then stays more or less constant 


after 2 : ~ 2. Total of ~ 4 x 10^ sources have formed 
in the simulation box of (75/i“^Mpc)^. Middle panel 
shows the resulting mass averaged densities of gas ki¬ 
netic energy, gas thermal energy and cosmic ray energy, 
which reflect mostly the averaged energy densities in¬ 
side clusters with high mass density. In clusters the 
cosmic ray energy dominated over the gas thermal en¬ 
ergy before z ~ 2 in our model, during which the gas 
thermal energy was still small. After z ~ 2, the ratio 
of cosmic ray to gas thermal energies has decreased. 
The cosmic ray energy density has become and stays 
~ 50% of the gas thermal energy density after z ~ I. 
With such amount of cosmic ray energy, the gas den¬ 
sity perturbation has reduced by ~ 20% at the present 
epoch, while the density perturbation of dark matter 
has been hardly affected, as shown in top panel. 

Fig. 4 shows the power of density perturbation at 
three different epochs. It has been decreased by ~ 25% 
at z = 2, by ~ 35% at z = 1, and by ~ 45% at z = 0 
on the cluster scale of ~ lh~^ Mpc, due to the dy¬ 
namical effects of cosmic ray pressure. However, the 
structures of scales larger than the cluster scale have 
been less disturbed, since sources form mostly at the 
highest density peaks within clusters. 

To study the dynamical effects of cosmic ray pres¬ 
sure on individual clusters, we have randomly identified 
several clusters in the computational box, and com¬ 
pared their properties in the cases with/without cosmic 
ray dynamics. Cosmic ray pressure provides additional 
support against infalling flows. As a result, the intra¬ 
cluster medium becomes less concentrated with cosmic 
rays, and the virial temperature reduces to lower val¬ 
ues. The clusters with Eqr ~ l/2Eth have been signif¬ 
icantly modified. Their X-ray bremsstrahlung emission 
and temperature decreased by ~ 50%. 

Further details of the simulations in this section will 
be reported in Ryu et al. (2002a). 

This work was supported by grant No. R01-I999- 
000-00023-0 from the Korea Science & Engineering 
Foundation. Simulations have been made through a 
support by the Grand Challenge Program of KISTI Su¬ 
percomputing Center. 
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Fig. 3. — Top: Temporal evolution of the root-mean-square density perturbations of dark matter and gas in the 512® 
ACDM simulations with (solid line) or without (dotted line) cosmic ray energy ejection from AGN-like sources. Middle: 
Temporal evolution of the mass-averaged energy per unit comoving volume. Ekin is the gas kinetic energy, Eth is the gas 
thermal energy, and Ecr is the cosmic ray energy. Bottom: Formation history of cosmic ray sources. Sources are assumed 
to form at 40 logarithmically spaced epochs after 2 = 10. 
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